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Abstract 

The exclusive production of one /i^/i^ pair in collisions of two ultra-relativistic nuclei is 
considered. We present a simple method for the calculation of the Born cross section for this 
process based on an improved equivalent photon approximation. We find that the Coulomb 
corrections to this cross section (corresponding to multi-photon exchange of the produced //^ 
with the nuclei) are small while the unitarity corrections (corresponding to the exchange of 
light-by-light blocks between nuclei) are large. This is in sharp contrast to the exclusive e'^e~ 
pair production where the Coulomb corrections to the Born cross section are large while the 
unitarity corrections are small. We calculate also the cross section for the production of one 
pair and several e'^e" pairs in the leading logarithmic approximation. Using this cross 
section we find that the inclusive production of /^^^^ pair coincides in this approximation with 
its Born value. 

PACS numbers: 
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I. INTRODUCTION 



Lepton pair production in ultra-relativistic nuclear collisions were discussed in numer- 
ous papers (see |l| for a review and references therein). For the RHIC and LHC colliders 
the charge numbers of nuclei = = Z and their Lorentz factors 71 = 72 = 7 are 
given in Table HI 



TABLE I: Colliders and cross sections for the lepton pair production 



Collider 


Z 


7 






RHIC, Au-Au 


79 


108 


36.0 


0.23 


LHC, Pb-Pb 


82 


3000 


227 


2.6 


LHC, Ar-Ar 


18 


3400 


0.554 


0.0082 



The cross section of one e^e^ pair production in Born approximation, described by 
the Feynman diagram of Fig. 1, was obtained many years ago P|. Since the Born cross 
section o^^^^ is huge (see Table Hj), the e^e~ pair production can be a serious background 
for many experiments. It is also an important issue for the beam lifetime and luminosity 
of these colliders 0. It means that various corrections to the Born cross section, as well 
as, the cross section for n-pair production, are of great importance. At present, there 
are a lot of controversial and incorrect statements in papers devoted to this subject. The 
corresponding references and critical remarks can be found in Refs. 0,0,0]. 

Since the parameter Za may be not small [Za ~ 0.6 for Au-Au and Pb-Pb collisions), 
the whole series in Za has to be summed in order to obtain the cross section with sufficient 
accuracy. The exact cross section for one pair production a\ can be split into the form 

= CTBorn + CTCoul + CTunit , (1) 

where two different types of corrections need to be distinguished. The Coulomb correction 
o"coui corresponds to multi-photon exchange of the produced with the nuclei (Fig. |2I); 
it was calculated in jj]. The unitarity correction cTunit corresponds to the exchange of 
light-by-light blocks between nuclei (Fig. E)); it was calculated in It was found in Q] 



p 



-p. 



FIG. 1: The Feynman diagram for the lepton pair production in the Born approximation 



FIG. 2: The Feynman diargam for the Coulomb correction 



and M that the Coulomb corrections are large while the unitarity corrections are small 
(see Table ITT|l . The results of p| were confirmed recently in 'f?] by a direct summation of 
the Feynman diagrams. 



FIG. 3: The Feynman diagram for the unitarity correction 



In this paper we present our calculations related to the exclusive and inclusive muon 
pair production. This process may be easier to observe experimentally than e~^e~ pair 
production described above. It should be stressed that the calculation scheme, as well as, 
the final results for the yU^/i^ pair production are quite different than those for the e+e~ 
pair production. 

In the next section we calculate the Born cross section for one yU^yU" pair production 
using the improved equivalent photon approximation with an accuracy about 5 %. In 
Sect. 3 we present the Coulomb and unitarity corrections to the exclusive production 



TABLE II: Coulomb and unitarity corrections to the e^e pair production 



Collider 


fCoul 
CBorn 


Cunit 
CBorn 


RHIC, Au-Au 


-25% 


-4.1% 


LHC, Pb-Pb 


-14% 


-3.3% 


LHC, Ar-Ar 


-1.06% 


-0.025% 



of one muon pair. In this section we also obtain the probabihty for the production of 
one muon and several electron pairs in collisions of nuclei at a given impact parameter p. 
This result allows us to calculate the inclusive production of one muon pair in Sect. 4. In 
Sect. 5 we calculate the cross section for two muon pair production in lowest order. The 
conclusion of the paper is given in Sect. 6. In the Appendix we consider the relatively 
simple process of one muon pair production by a real photon off a nucleus Zp^ pT 

which serves as a good test of the approach used. A short preliminary version of this 
paper was published in 

II. BORN CROSS SECTION FOR ONE fi+ fi' PAIR PRODUCTION 

The production of one pair 

+ ^ + + (2) 

in the Born approximation is described by the Feynman diagram of Fig. 1. When two 
nuclei with charges Zic and and 4-momenta Pi and P2 collide with each other, they 
emit equivalent (virtual) photons with the 4-momenta gi, g2, energies cji, UJ2 and their 
virtualities Q\ = —qf, Ql = — g|. Upon fusion, these photons produce a fi^p^ pair with 
the total four-momentum gi + ^2 and the invariant mass squared W'^ = {qi + q2Y- Besides 
this we denote 

{Pi + P2f = AE^ = AM"^ -f^ , a ^1/137 
and use the system of units in which c = 1 and h = 1. 
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A. General formulae 



The Born cross section of the process (j21) can be calculated to a good accuracy using the 
equivalent photon approximation (EPA) in the improved variant presented, for example, 
in Ref. ^. Let the numbers of equivalent photons be dni and dn2- The most important 
contribution to the production cross section stems from photons with very small virtuali- 
ties <^ /i^ where fi is the muon mass. Therefore to a good approximation, the photons 
move in opposite directions, and W"^ ~ 4ci;iCo'2. In this very region the Born differential 
cross section da-^ for the process considered is related to the cross section cr^^ for the 
process with real photons: 77 fi^^^ by the equation 

daB = dnidn2 d(j^^{W^) , W'^ ^ 4ujiuj2 ■ (3) 



The number of equivalent photons are (see Eq. (DA) in Ref. 



where 



UJ. 



2 



Qi > Qi mill — 2 

7 

and F{Q'^) is the electromagnetic form factor of the nucleus. It is important that the 
integral over converges rapidly for > l/i?^, where 

R= 1.2A^/3 fm (6) 

is the radius of the nucleus with A ^ M/ irip the number of nucleons (i? ~ 7 fm, ^ 28 
MeV for Au and Pb). Since Q^in ~ 1/-^^) the main contribution to the cross section is 
given by virtual photons with energies 

u^i<l/R- (7) 

Therefore, we can use the spectrum of equivalent photons neglecting terms proportional 
to Ui/Ei given by: 

After the transformation 

dijj\ dijj2 dui dW"^ 
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we cast the cross section in the form 



7r2 UJi V l^QlJ Ql 



where a;2 ~ W^/IAlvi). 



B. Leading logarithmic approximation (LLA) 

Before using the calculation scheme above, it is instructive to present a rougher but 
simpler approximation — the so called leading logarithmic approximation (LLA). In the 
LLA, the equivalent photon spectrum as a function of photon energy dni{uji) is obtained 
after integrating dni{u!i, Ql) over Ql in the region between 

QLin < Ql < ^IR^ (11) 

which leads to 

dmiui) ^ In . 12 

The restriction mm ^1/-^^ corresponds to the integration interval 

a = ^<a,,<i,= l (13) 

which gives 

Since a^^{W^) ^ (Ana'^/W^) In (W'^/iJ,^) for large values of 1^ » the main contribution 
to the Born cross section comes from the region of small values of W near the threshold. 
Therefore, within logarithmic accuracy wc replace W by some fixed value Wq ~ in the 
lower bound a. After that the integral over gives 

and further integration over cui leads to 

-^ln—ln^ = -L\ 16 
uji u!{ 3 
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where 



As a result, we obtain 



LLA 



28 {ZiaZ2af 



277r /i^ . . 

in accordance with the result of Landau & Lifshitz The accuracy of the LLA depends 
on the choice of the value for Wq. If we use for numerical estimations Wq = 3/i, then the 
accuracy of the LLA for the colliders discussed is about 15%. 

The same result can be obtained in the framework of the impact parameter dependent 
representation, which will also be useful later. For this aim we introduce the probability 
for muon pair production -Pb(p) in the collision of two nuclei at a fixed impact parameter 
p. For 7 ^ 1, it is possible to treat the nuclei as sources of external fields and to calculate 
Pb(p) analytically using the same approach as in Ref. Id]. The Born cross section as can 
then be obtained by the integration of -Pb(p) over the impact parameter: 

aB= [ PB{p)d^p. (19) 



We calculate this probability in the LLA using Eq. ^ with 

dni = ^, cu, i?<p, (20 

TT^ OJi Pi R Ui 

where Pi is the impact parameter of z-th equivalent photon with respect to the i-th nucleus. 
This allows us to write the above probability in the form 

Pb(p) = j dn.dn, 6ip, -p^-p) a,,{W^) = ^ ^-^J^^ ^(.P) ■ (21) 
Depending on the value of p two different forms for $(p) need to be used: 

<|.(p)= (41n^ + ln^) In^ for i? « p < ^ , (22) 

(2 \ 2 2 
In^L ) for 2<p<^. (23) 

p'pRj p-^ p^R ^ ^ 

Note that the function $(p) is continuous at p = '-//p together with its first derivative. 

As expected the integration of -Pb(p) over p in the region R < p < 7^/(/i^/2) gives back 

the result in (fTHj) . 

To prove ()2Hl - (f^ . we make the transformation given in together with the inte- 
gration over W"^ according to (jl5|) . This gives 

„ , , 14 (ZiaZ2ay P'^ dooi f d^pi J ^ \ /7CJ1 , ,V^^, 
Pb{p) = —, / ' / — / 2r --piUM-^- P-Pi 24 



where "^{x) is the step function. The main contribution to this integral is given by two 
regions: R <^ pi ^ p and R ^ \p — pi\ = p2 ^ p- Moreover, the two regions in the oji 
integration with p < ui < ■j/ R and p'^R/j < Ui < p give the same contributions. As a 
result we get 

= 2(Ji + J2) , (25) 



— [ —^{'^/uJi-pi)^{-iUJi/p^-p) 



Ji 

^1 Jr Pi 

J2 = r""— r^^(7M-p)^(7^i//.^-P2). 

Jt, Jr P2 

Next we consider the two regions of p. 

In the region of relatively small impact parameters, R <ti p < 'j/ p, the second step 
function does not impose any limitations, therefore. 



^ -"'d^ I'dp^^ r'-,j^ r'-*l = ,„X,„^ + i(,„iL 

^1 Jr Pi J-y/p Jr Pi pp R 2 \ R 
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^1 Jr Pi pp R 
Summing up, we obtain (f^^ 

In the region of relatively large impact parameters, 'y/p<P'^ 7^/(/i^-R), we have 

^ 7^.^/7 ^1 Jr Pi 2 V p^Rp) ' ^ 



therefore, the sum gives 



culations based on the exact 
. We find good agreement 



We compare Eqs. (j?T|) - ()23|) for with the numerical ca 
matrix element calculated with the approach as outlined in 
for Pb-Pb collisions: the discrepancy is less then 10 % at pp > 10 and it is less then 15 % 
at pp > 2pR = 7.55. 

C. More refined calculation 

In the calculation below we use for the form factor of the nucleus the simple approxima- 
tion of a monopole form factor, which corresponds to an exponentially decreasing charge 



distribution, whose mean squared radius \/Jr^ is adjusted to the experimental value: 



For lead and gold, the parameter is A 80 MeV. This approximate form of the form 
factor enables us to perform some calculations analytically, which otherwise could only 
be done numerically. 

The equivalent photon spectrum dniiuji) is obtained after integrating driiiuji, Qf) over 
(the upper limit of this integration can be set to be equal to infinity in a good approx- 
imation, due to the fast convergence of the integral at > A^): 

dnM)-^f(^)'-^. (27) 



Here the function 



/(x) = (l + 2x^)ln(^ + l) -2 (2^ 



is large for small values of x, 



/(x) ^In^ - 2 = ln-^ at x< 1, (29) 

ex r 



but drops very quickly for large x in accordance with Eq. ((7j): 



/(x)<^ for x>l. (30) 



Finally we obtain 



where 



^max 



(32) 



Since uOi < E and UJ1UJ2 ~ we have cjmin ~ ^J''^ / E and cUmax = E. However, due to the 
fast decrease of /(x) for x > 1 one can extend these limits up to dJmin = and tUmax = 00 
without any lack of accuracy, therefore, 

00 

G{W') = 2 I fix,) /(X2) dy , xi,2 = ^ (33) 


A numerical evaluation of the integrals in Eqs. (jHH), (jH^ yields the function J{'jA/fj,) 
presented in Fig. |3] 

Corrections to the photon spectrum are represented by terms in Eqs. (jSl), (0) of the 
order of Q^/W^ (see Eqs. (E.l) in Ref. Q]), which are dropped before the integration over 
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FIG. 4: The function J{-fA/fi) from Eq. ^ 



Qj is done. After the integration with weight 1/Qf the relative value of these corrections 
becomes of the order of 



(34) 



Thus, for the collisions considered here one can estimate the accuracy of the calculations 
on the level rji ~ 5%. Another test of accuracy of the approach used is given in Appendix. 
Note that in the LLA the function GiW"^) is just 



In 



n 3 



eiy/2A)2 



In 



{WR/2y 



(35) 



in accordance with Eq. p6|) (taking into account that A/e ~ ^/R)- 



III. COULOMB AND UNITARITY CORRECTIONS 

The Coulomb correction corresponds to Feynman diagram of Fig. 2. Due to the re- 
striction of the transverse momenta of additionally exchanged photons to the range below 
l/R, the effective parameter of the perturbation series is not (Za)^ but (Za)^ /{Rfi)^. 
In addition, the contribution of the additional photons is suppressed by a logarithmic 
factor. Indeed, the cross section for two-photon production mechanism is proportional to 
L^, while the cross section for the multiple-photon production mechanism is proportional 
only to L^. Therefore, the real parameter describing the suppression of the Coulomb 
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correction is of the order of 

which corresponds to Coulomb corrections of less then 1 %. The example considered in 
Appendix confirms this estimate. 

The unitarity correction (Xunit to the one muon pair production corresponds to the 
exchange of light-by-light blocks between the two nuclei (Fig. 3). We start with a more 
general process — the production of one fi~^fi~ pair and n electron-positron pairs (n > 0) 
in a collision of two ultra-relativistic nuclei 

Zi + Zi + Z2 + i^+^i- +n{e+e~) (37) 

taking into account the unitarity corrections, which correspond to the exchange of the 
blocks of light-by- light scattering via the virtual lepton loops. The corresponding cross 
section dai^n can be calculated by a simple generalization of the results obtained in 
for the n-pair process without muon pair production: Zi -\- Z2 —>■ Zi + Z2 + n(e~^e~). 
This multiple pair production process was studied in a number of papers, see [l[ for a 
review. It was found that the probability is to a good approximation given by a Poisson 
distribution with the deviations found to be small. Indeed, it is not difficult to show that 
the ba* equation, fo. tKe .atte. p.oeess should be .od,«ed a. follows. In E, (28) of Q 
the additional factor 

5^(p,r„+i)e-^[^"('')/2+iv,-(p)] (38) 

appears under the integral, where L = In (7172) and the functions B'^, A'^ and <y9^ are the 
same as the functions B, A and if in Eq. (27) of but with the replacement of electrons 
by muons. A. a .suit. E., ,31) of . .placed by 

^ = L^;(p)L?dlP!e-"f <">-"■("', (39) 

where LA'^{p) ^ -Pb(p) is the probability for one muon pair production in the Born 
approximation, as discussed in Sect. 2.2. In the region of interest, p > 2R, the function 
Ai{p) is small, 

L<(p)«l, A5'(p)«A(p), (40) 
therefore, we can rewrite ()39|1 in the simpler form 

^ = P,Up) , PiUp) = PM e-^('') , (41) 
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where ne(p) = LAi{p) is the average number of e'^e~ pairs produced in colhsions of the 
two nuclei at a given impact parameter p. The result that the probabilities for the different 
processes factories is due to the independence of the individual processes. For a general 
discussion of the validity of this factorization together with possible violations we refer to 

0. 

In particular, we get the cross section for the exclusive one pair production 

including the unitarity correction as 

^1+0 = I PB(p)e-^^(^)f/V- (42) 
This expression can be rewritten in the form 

C^l+O = 0"B + O^unit , (43) 

where 

crB = j Pb{p) d^P (44) 
is the Born cross section discussed in Sect. 2 and 

cTunit = -/[!- e"'''^^''^ Mp) d'p (45) 

corresponds to the unitarity correction for the one muon pair production. 

A rough estimation of (Xunit can be done as follows. The main contribution to cTnmt 
comes from the region 

i? < p < l/rrie (46) 

in which the function ne(p) ~ ?^e(2-R) and the integral ()45|) can be calculated in LLA. It 
gives 

f^unit ~ - — [1 - e ^ Junit , (47) 

where 

Junit = 6/ $(p)^. (48) 



hR P 

As a result, we find Cumt ~ —1.2 barn for the Pb-Pb collisions at LHC, which corresponds 
approximately to (—50)% of the Born cross section. 

It is seen that unitarity corrections are large, in other words, the exclusive production 
of one muon pair differs considerable from its Born value. 
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IV. INCLUSIVE PRODUCTION OF ONE fi+fi- PAIR 



The experimental study of the exclusive muon pair production seems to be a very 
difficult task. Indeed, this process requires that the muon pair should be registered 
without any electron-positron pair production including emitted at very small angles. 
Otherwise, the corresponding cross section will be close to the Born cross section. 

To prove this, let us consider the process ()37|) . whose probability is given by Eq. (PT|) . 
The corresponding cross section is 

ai+n = J Pi+n{p)d?p. (49) 

It is clearly seen from this equation that after summing up over all possible electron pairs 
we obtain the Born cross section 

oo 

CTl+n = 0"B • (50) 

n=0 

Therefore, there is a very definite prediction: the inclusive muon pair production coincides 
with the Born limit. This direct consequence of calculations, which take into account 
strong field effects, may be easier to test experimentally than the prediction for cross 
sections of several e^e~ pair production. 

V. TWO MUON PAIR PRODUCTION 

The cross section of the process 

Z1Z2 Z1Z2 + ii^p~ii^ix~ (51) 

can be calculated in lowest order in a according to 

^2 = ^1 [PB{p)fd^P (52) 

with the integration region p > 2R. But in this region the probability Pb{p) is given to a 
good accuracy by Eqs. (j^Hjl - lj^^j) . From this we get a2 = 1.24 mbarn for Pb-Pb collisions 
at LHC. 

VI. CONCLUSION 

The exclusive production of one fi~^fi~ pair in collisions of two ultra-relativistic nuclei 
is considered. We present a simple method for the calculation of the Born cross section 
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for this process. 

We found that the Coulomb corrections to this cross section are small, while the unitar- 
ity corrections are large. This is in sharp contrast to the exclusive e~^e~ pair production 
where the Coulomb corrections to the Born cross section are large while the unitarity 
corrections are small. 

We calculate also the cross section for the production of one ^~^fi~ pair and several 
e~^e~ pairs in LLA. Using this cross section we found that the inclusive production of 
fj^~^fi~ pair coincides in this approximation with its Born value. 

Let us discuss the relation of the cross sections obtained for the muon pair production 
with the the differential cross section of the e'^e" pair production in the region of large 
transverse momenta for the e^, for example at p±± > 100 MeV. It is clear that for the 
e~^e~ pair production in this region, the situation is similar to the case considered for 
yU+yU" pair production. 

We expect that the inclusive production of a single e^e~ pair with large transverse mo- 
menta of (together with several unobserved e~^e~ pairs in the region of small transverse 
momenta of of the order of me) coincides with the Born limit. 



To tests the approach used in Sect. 2.3, we consider the simpler case of the muon pair 
production by a real photon with the energy u off an nucleus 



This cross section was calculated by Ivanov and Melnikov in Ref. |11| using the same 
expression for the form factor of the nucleus and assuming A^/(2/i)^ <^ 1. The 
corresponding formula for the Born contribution and the first Coulomb correction is 



Appendix 



'yZ Zn II 



(53) 




(54) 



where 




(55) 



Therefore, the relative magnitude of the Coulomb correction is given by 



C2 



(56) 
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which confirms the estimate in ()36|) . 

In the equivalent photon approximation, the cross section is given by 

da^^^ = dn2(T^^{W^) (57) 

which has the form 

EPA Z^a p dW^ ^ fW^\ ,,,.2^ 



The main contribution to this integral is given by the region near the lower limit, where 
the argument of the function / is small and therefore / can be replaced by its approximate 
expression ^2^ : 

/ ^ =2 1n^-2-2 1n^. (59) 
After that the cross section can be calculated without difficulties as 



^EPA o 



IT 

where / is given by (jl5p and 



28 Z^a^ 



9 fx^ 



I , (60) 



f^dW^ .X , (43- 281n2)7ra2 

Comparing this expression with the one of ()54|) . we find that those terms, which are 
omitted in the EPA, have a relative magnitude of the order of 

^1 = y , (62) 
this expression confirms the estimate (j34p . 
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